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We describe a formalism to calculate form factor and charge density distribution of the 
pion in the chiral limit using the holographic dual model of QCD with hard-wall cutoff. 
We introduce two conjugate pion wave functions and present analytic expressions for these 
functions and for the pion form factor. They allow to relate such observables as the pion 
decay constant and the pion charge electric radius to the values of chiral condensate and 
hard- wall cutoff scale. The evolution of the pion form factor to large values of the momentum 
transfer is discussed, and results are compared to existing experimental data. 
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INTRODUCTION 



During the last few years applications of gauge/gravity duality [l[ to hadronic physics attracted 



a lot of attention, and various h olog raphic dua 



(see, e.g., 



11 



12 



models of QCD we re p roposed in the literature 



14, 15 



,13 , £ Q, Q, I21Q. These models 



were able to incorporate such essential properties of QCD as confinement and chiral symmetry 
breaking, and also to reproduce many of the static hadronic observables (decay constants, masses), 
with values rather close to the experimental ones. Amongst the dual models, a special class is 
the so-called "bottom-up" approaches (see, e.g., Q, Q, 0, []]), the goal of which is to reproduce 
known properties of QCD by choosing an appropriate theory in the 5-dimensional (5D) AdS bulk. 
Within the framework of the AdS/QCD models, by modifying the theory in the bulk one may try 
to explain/fit experimental results in different sectors of QCD. 



In the present paper, we will be interested in the hard- wall AdS/QCD model Q, Q] , where 
the confinement is modeled by sharp cutting off the AdS space along the extra fifth dimension 
at a wall located at some finite distance z = Zq. In the framework of this hard- wall model, it is 
possible to find form factors and wave functions of vector mesons (see, e.g., [23]). To reproduce 
the general features of the spectrum for the higher states ("linear confinement"), a soft- wall model 



ti {{J. 



24]. 
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In general, the vector sector is less sensitive to the infrared (IR) effects, since this symmetry is 
not broken in QCD. However, the axial-vector sector appears to be very sensitive to the particular 
way the chiral symmetry is broken or, in other words, to the bulk content and the shape of the IR 
wall 0]. 

In this respect, one of the interesting objects to study in the holographic dual models of QCD 
is the pion. The properties of the pion were studied in various holographic approaches, (see e.g. 



Refs. [a, 
papers [if 




221]). In particular, the approach of Ref. [a] (see also recent 



211 ]) managed to reproduce the (Gell-Mann-Oakes-Renner) relation 



m 



between the quark mass m q and mass of the pion m n and also the g pn7T coupling (the coupling 
between p meson and two pions). In Ref. the solution of the pion wave- function equation was 
explicitly written for the m„ = limit. 

n 

In this paper, working in the framework of the model proposed in [6( (hard- wall model), we 
describe a formalism to calculate the form factor and wave functions (and also the density function) 
of the pion. Since the fits of Ref. [fi] give a very small m q ~ 2MeV value for the explicit chiral 
symmetry breaking parameter m q , we consider only the chiral limit m q = of the hard- wall 
holographic dual model of two-flavor QCD. Resorting to the chiral limit allows us to utilize one 
of the main advantages of AdS/QCD - the possibility to work with explicit analytic solutions of 
the basic equations of motion. Expressing the pion form factor in terms of these solutions, we are 
able, in particular, to extract and analyze the behavior of the pion electric radius in various regions 
of the holographic parameters space. On the numerical side, we come to the conclusion that the 
radius of the pion is smaller than what is known from experiment. However, we suggest that, as 
in case of the radius of the p meson, smoothing the IR wall may increase the pion radius. 

In our analysis, we introduce and systematically use two types of holographic wave functions 
<&(z) and ^f(z), which are conjugate to each other and basically similar to the analogous objects 



introduced in our papers 



23 



241 ] . where we studied vector mesons. 



The paper is organized in the following way. We start with recalling, in Section II, the basics of 



the hard- wall model and some results obtained in Ref. [y], in particular, the form of the relevant 
action, the eigenvalue equations for bound states and their solutions. In Section III, we describe 
a formalism for calculating the pion form factor and express it in terms of the two wave functions 
mentioned above. In Section IV, we discuss the relation of our AdS/QCD results to experimental 
data. We express the values of the pion decay constant and the pion charge radius in terms of 
the fundamental parameters of the theory and study their behavior in different regions of the 
parametric space. At the end, we study the behavior of the pion form factor at large momentum 



3 



transfer. Finally, we summarize the paper. 



II. PRELIMINARIES 



In the holographic model of hadrons, QCD resonances correspond to Kaluza-Klein (KK) exci- 
tations in the sliced AdSs background. In particular, vector mesons correspond to the KK modes 
of transverse vector gauge field in this background. Since the gauge symmetry in the vector sector 
of the H-model is not broken, the longitudinal component of the vector gauge field is unphysical, 
and only transverse components correspond to physical mesons. Similarly, the axial-vector mesons 
are the modes of the transverse part of the axial- vector gauge field. However, because the axial- 
vector gauge symmetry is broken in the 5D background, the longitudinal components have physical 
meaning and are related to the pion field. This should be taken into account if we want to treat 
the pion in a consistent way. 



A. Action and Equations of Motion 

The standard prescription of the holographic model is that there is a correspondence between 
the 4D vector and axial-vector currents and the corresponding 5D gauge fields: 

J^(x) = q(x)^t a q(x) - V£(x, z) (1) 
J%»(x) = q{xh^ 5 t a q(x) A*(x,z) , 



where t a = cr a /2, (a = 1, 2, 3 and <r a are usual Pauli matrices). 

In general, one can write A = A± + An, where A± and An are transverse and longitudinal 
components of the axial-vector field. The spontaneous symmetry breaking causes Au to be physical 
and associated with the Goldstone boson, pion in this case. The longitudinal component may be 
written in the form: ^.^-..(x, z) = dMip a (x, z). Then ip a (x, z) corresponds to the pion field. Physics 
of the axial-vector and pseudoscalar sectors is described by the action 



S£ dS = Tr f d 4 x dz \-AD M X)\D M X) + lx*X - -\-A MN A M 



(2) 



where DX = dX — iA^X + iXA^, (Aum = V ± A) and X(x,z) = v(z)U(x, z)/2 is taken as a 
product of the chiral field U(x, z) = exp (2it a 7r a (x, z)) and the function v(z) = m g z+az 3 containing 
the chiral symmetry breaking parameters m q and a, with m q playing the role of the quark mass 
and a that of the quark condensate. Expanding U(x, z) in powers of 7r a gives the relevant piece of 
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the action 



S 



A (2) 
AdS 



Tr / <Fx dz 



1 aMN a i V 2 ( z ) f A g r, a\2 

A Amn + „ 5 \ A M - ) 



±9h 



2z 3 



(3) 



This Higgs-like mechanism breaks the axial- vector gauge symmetry by bringing a z-dependent mass 
term in the ^4-part of the lagrangian. Varying the action with respect to the transverse part of the 
axial-vector gauge field A°[_ (x, z) and representing the Fourier image of A^^x, z) as A\_ (p, z) we 
will get the following equation of motion 

1 



-8 A a 



2„,2 Aa 



o , 



(4) 



that determines physics of the axial- vector mesons, like A±. The axial- vector bulk-to-boundary 
propagator A(p, z) is introduced by the relation Aj_ (p, z) = A(p, z)A®(p). It satisfies Eq. (j3|) with 
boundary conditions (b.c.) A(p, 0) = 1 and A'(p,zq) = 0. Similarly, variation with respect to the 
longitudinal component d^ip a gives 



z 3 d z 



i 



d z r - 55 v(^-7r a ) = o. 



(5) 



Finally, varying with respect to A z produces 

P 2 z 2 d z r - g!v 2 d z 7T a ■- 
) and 



0. (6) 
} with b.c. d z ifj(zo) = 0, -0(e) = and 



The pion wave function is determined from Eqs. 
vr(e) = 0. 

Within the framework of the model of Ref. (|, it is possible to derive the Gell-Mann-Oakes- 
Renner relation m 2 ~ m q producing massless pion in the m q = limit. Taking p 2 = m 2 in Eq. ([6 
gives 



8 Z TT 



(7) 



A perturbative solution in the form of m 2 expansion was proposed in Ref [g], with ip{z) = .4.(0, z) — l 
in the lowest order. Then it was shown that, in the m q — > limit, tt(z) tends to — 8(z — zq) or, 
roughly speaking, tt = —1 in this limit. Since our goal is to calculate the pion form factor in the 
chiral limit, this approximation will be sufficient for us. 



B. Two-Point Function 



The spectrum in the axial-current channel consists of the pseudoscalar pion (0| J^\tt(p)) = if-KP a 
and axial-vector mesons (0\J^\A n (p,s)) = F^ n e"(p, s), where F^ n correspond to the n th axial- 
vector meson decay constant (and we ignored the flavor indexes). Thus, the two-point function for 
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the axial-vector currents has the form: 



where the meson polarization tensor is given by 
The representation for the two-point function can be also written as 

( Ja(p)4(-p) ) = pV4 + f-^ a/3 + ^) E 2 + ( non p° le terms )> ( 10 ) 

P V P / „ P - M A,n 

in which the second term on the rhs is explicitly transverse to p. 

As noted in Ref. fjj], using holographic correspondence one can relate the two-point function to 
[d z A(p, z)/z] z= q and derive that 

fl = -^(-d z A^ z )) . (11) 

For large spacelike p 2 , Eq. @ gives the same solution as in case of vector mesons, and the same 
asymptotic logarithmic behavior, just as expected from QCD. 

C. Pion Wave Functions 

The longitudinal component of the axial- vector gauge field was defined as Am = dip. In the 
chiral limit, when p 2 = m 2 = 0, we have d z ir = 0, and the basic equation for ip, Eq. (J5j) can be 
rewritten as the equation 

z 3 d z d z *j - g 2 5 v 2 * = (12) 

for the function = ip — ir. In the chiral limit, when tt(z) — ► — 1, the value of ^(e) tends to 1 
as e — ► 0. This value and the b.c. ^'(zo) = are the same as those for A(p, z) and, furthermore, 
Eq. (fT2|) coincides with the p 2 = version of equation for A(p, z). Hence, the solution for ^(z) 
coincides with A(0, z): 



and we may write 



*(z)=A(0,z), (13) 



35 \ z / z=e^0 



In our analysis of p-meson wave functions in Refs. [23|, [2J], we emphasized that it makes sense 
to consider also the conjugate functions $(z) ~ *$>'(z)/z of the corresponding Sturm-Liouville 
equation. As we observed, they are closer in their structure to the usual quantum mechanical 
bound state wave functions than the solutions of the original equation. In the pion case, it is 
convenient to define the <I> function as 

•w = -iBf (I5) 



It vanishes at the IR boundary z = zq and, according to Eq. (jTTj) . is normalized as 

$(0) = 1 (16) 
at the origin. Note also that using Eq. (fT2|) we can express * as derivative of <I>: 

*{z) = -£^-d t ${ k z) . (17) 

III. EXTRACTING PION ELECTROMAGNETIC FORM FACTOR 

A. Three-point function 

To obtain the pion form factor, we need to consider three-point correlation functions. The 
correlator should include the external EM curr ent J^(0) and currents having nonzero projection 
onto the pion states, e.g. the axial currents Jg a (xi), J§l(x2) 

% a p(pi,p 2 ) = J d 4 x, J d A x 2 e'w-**** (0\Tjl p (x 2 )j; l (0)J 5a (xiM , (18) 

where pi,p 2 are the corresponding momenta, with the momentum transfer carried by the EM 
source being q = P2 — Pi (as usual, we denote q 2 = —Q 2 , Q 2 > 0). The spectral representation for 
the three-point function is a two-dimensional generalization of Eq. ([8]) 

f 2 F (Q 2 ) 

T^ a P(pi,p2) = PiP^ipi + P2Y— — 1" (transverse terms) + (nonpole terms), (19) 

P1P2 n,m 

where the first term, longitudinal both with respect to pf and p^ contains the pion electromagnetic 
form factor F 7T (Q 2 ) 

Mpi)|J?(0)br(p2)> = F n (q 2 ){p 1+ p2)^ , (20) 

(normalized by F n (0) = 1), while other pole terms contain the contributions involving axial- vector 
mesons and are transverse either with respect to pf or p^, or both. Hence, the pion form factor 
can be extracted from the three-point function using 

PlaP2(3T^( Pl ,p2)\ p 2 =0tpl=0 = {pi+P2YflF^Q 2 ) . (21) 
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B. Trilinear Terms in F 2 Part of Action 

To obtain form factor from the holographic model, we need the action of the third order in the 
fields. There are two types of terms contributing to the pion electromagnetic form factor: |DX| 2 
term and F 2 terms. Let us consider first the contribution from F 2 terms. They contain VVV, 
VAA and AVA interactions and may be written as 

^Adsls = ^ Tr / d\ dz - {V MN [V M , V n ] + V MN [A M , A N ] + A MN [V M , A N ]) , (22) 

where Vmn = QmVn - 9nVm and Amn = 9mAn — djvA/Vf- Taking V z = A z = gauge, we pick 
out the part of the action which is contributing to the 3-point function {J§ a J^J^n)'- 

W 3 = ^Tr [ d 4 x dz - (V^[A», A v ] + A^[V», A v \) . (23) 

Introducing Fourier transforms of fields, we define, as usual, VAq, z) = V^(q)V(q, z) for the vector 
field, where V^{q) is the Fourier transform of the 4-dimensional field V^x) and V(q, z) is the 
bulk-to-boundary propagator satisfying the equation 

zd z (^d z V(q,z)\ +q 2 V(q,z) =0 (24) 

with b.c. V(q, 0) = 1 and d z V(q, Zq) = 0. It can be written as the sum 

Vh.*)=»f;^^r <*> 

m=l 



involving all the bound states in the g-channel, with M m being the mass of the mth bound state 
and ipm(z) its wave function given by a solution of the basic equation of motion in the vector sector. 

The projection (|21|) picks out only the longitudinal part A^(p, z) of the axial- vector field. 
Taking into account that AW (x, z) = d fl ip(x, z), we may write 

A^(p,z)=ip^ a (p,z) , (26) 

where A^(p,z) and tp a (p,z) are the Fourier transforms of A$(x,z) and vp(x,z), respectively. 
Furthermore, there is only one particle in the expansion over bound state in this case, namely, the 
massless pion. Thus, we have A% (p, z) = Ajf {p)tp(z) and, therefore, 

4 a (p,z) = -K A t*<PW*) ■ (27) 
This allows us to rewrite AS (p, z) in the form 

A 1(P^) = ^ A 1(P)^) (28) 



8 



involving the longitudinal projector p a p^/p 2 and the pion wave function tp(z), which is the solution 
of the basic equation ([5]) . Using this representation and making Fourier transform of W3 gives 



(-abc 



d A u d A v d A w 



(27T) 



12 



i{2Tr) 4 5 w (u + v + w) 



U 2 V 2 



(29) 



A\ a {u)Al p {v) w,V v a {w)- Wu V°{w) / dz -V(w,z)^(z) 



1 



(notice that the second term in Eq. (|23p vanishes for longitudinal axial- vector fields). Varying this 
functional with respect to sources produces the following 3-point function: 



PIP2 



2si 



dz - V{q,z)4> 2 {z) 



(30) 



where, anticipating the limit p\ — > 0, p| — ► 0, we took (pig) = —(p2l) = —Q 2 /2 m the numerator 
factors. Now, representing (^ a (g)^f J 4 ife (Pi)^j 3 j4)C (-?'2)) = z(27r) 4 <5( 4 )(g + p x - p 2 ) ^abc^ al3 (pi,P2) 
and applying the projection suggested by Eq. (|2"T]). we will have 



lim lim pi a p2pT^( Pl ,p 2 ) = -L (pi +p 2 fQ 2 J(Q) , 
p?^o p^o 2^ 

where J(Q) is the dynamic factor given by the convolution 

" z ° dz 



J(Q) 



J(Q,z)^(z) . 



(31) 



(32) 



C. Dynamic Factor and Wave Functions 



The vector bulk-to-boundary propagator J{Q,z) = V(iQ,z) for spacelike momenta, entering 
into the dynamic factor J(Q), satisfies the equation 



zd z ( -d z J(Q,z) ) =Q z J(Q,z 



with b.c. i7(Q, 0) = 1 and d z J{Q, Zq) = 0. Its explicit form is given by 

K (Qzo)' 



J(Q,z)=Qz 



Kx(Qz) + h(Qz) 



One can easily see that S(0, z) = 1. Combining all the factors, we get 



fM F2) (Q 2 ) 



1 



2 Q 2 



dz 



J(Q,z)^ 2 (z) . 



igi Jo z 

Integrating by parts and using equations of motion both for J and tp gives 

2 ,2„,2 



1 



d\ fn J 



dzzJ{Q,z) 



d z tp 



+ 



^V(V'-vr) 



(33) 



(34) 



(35) 



(36) 
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We need also to add the Vtttt contribution from the |.DX| 2 term of the AdS action ([2]). It is 
generated by 

V(z) 



fi 1 AdS I ^ 7r7r — ^abc I d X dz 



iA a M -d M ir a )K b V c * 

&~ 

and its inclusion changes tp(ip — it) into {ip — it) 2 in Eq. (|36p . The total result (see also Ref. [(j) 
may be now conveniently expressed in terms of the = ip — tt wave function 



(37) 



FAQ 2 ) 



1 



dz zJ(Q,z) 



+ 9 ^\z) 



(38) 



Using equation of motion for ^(z), one can see that the expression in square brackets coincides 
with 

- z d z U{z)- z d z *{*)\ =-gtfl- z d z {*{z)<$>(z)) , 

and write the form factor as 

FAQ 2 ) = - r dz <?(Q> z ) d * (^( z ) ®( z )) • ( 39 ) 



This representation allows one to easily check the normalization 



fao) 



dz d z (v(z) = (0) $(0) = 1 



(40) 



where we took into account that J~(0,z) = 1 and &(zo) = 0. We can also represent our result for 
the pion form factor as 



FAQ' 



dz z J(Q, z) 



a 



9 2 5 & 2 (z) + ^z'^(z) 

Jit 



dz z J(Q, z) p(z) 



and interpret the function p(z) as the radial distribution density, as it was done in Refs. [2J; 



(41) 



Note that keeping only the first term in square brackets gives an expression similar to our result 



23j for the p-meson form factor 



Fn(Q 2 ) 



-0 



dzzJ(Q,z)\M*)\' 



(42) 



in terms of the function <j>\ conjugate to the solution of the basic equation of motion. The value 
of (f>i(z) at the origin is proportional to the p-meson decay constant f p /m p = g p (experimentally, 
g p xp « 207 MeV [251]). namely, 4>i(0) = g§g p . Thus, the pion wave function g^f v ^{z) = (p w (z) is a 
direct analog of the p-meson wave function (fti(z). Main difference is that, in the pion case, there 
is also the second term in the form factor expression. The latter, in fact, is necessary to secure 
correct normalization of the form factor at Q 2 = 0. In Eq. (|38p . this term is written in terms of 
the *$>(z) wave function, but using Eq. (|17p we can rewrite it also in terms of &(z) or (j)A z )'- 



p(z) = 4>l (z) + 



1 



9h 2 



1 



d z $A z \ 



(43) 
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IV. WAVE FUNCTIONS AND FORM FACTOR 
A. Structure of Pion Wave Functions 



Explicit form of the \& wave function follows from the solution of Eq. (jl2f) : 

(44) 



*(*) = *r[2/3] (|) 1/3 



t ( 3\ t ( 3\ ^2/3 ( a2; o) 

i_ 1/3 (az 3 ) - J 1/3 (az 3 ) v , 3 ' 

i_ 2 /3 [ az o) 



where a = w 1.481 <r (recall that 55 = \/2~7r, see e.g. Ref. 



24|). As a result, <&(z) is given by 



1 (\ \ / a 4\l/3 



, <jx ^2/3(«^) 



-/ 2/3 (az 3 ) + J_ 2/3 (az 3 ) 



J -2/3 M) 



(45) 

This formula, combined with Eq. (116f) . establishes the relation 



,2 o 9 i/3 £[2/3] J 2/3 (^o 3 ) « 2/3 , . 

/?r " r[l/3] J_ 2/3 (az 3 ) ^ (45) 

for /„- in terms of the condensate parameter a and the confinement radius Zq . Since a appears in 
the solutions only through a, we will use a in what follows. Note also that a 1 / 3 ~ 1.14a 1 / 3 . 

Realizing that the equations of motion for the vector sector in this holographic model are not 
affected by the chiral symmetry-breaking effects expressed through the function v(z), it is natural 
to set the value of zq from the vector sector spectrum, i.e., by the p-meson mass. The numerical 
value of zq (call it Zq) is then Zq ~ 1/323 MeV. As given by Eq. (|46p . f n looks like a rather 
complicated function of two scales, zq and a. Note, however, that the ratio l2/z{ a ) / I-2/z{ a ) is 
very close to 1 for a ^ 2 and practically indistinguishable from 1 for a ^ 3. Hence, for sufficiently 
large values of the confinement radius, zq > 1/a 1 / 3 , the value of f w is determined by the value of 
a alone. This limiting value of /„■ is given by 



/ttU-oo -1 g5 y r[1/3] - 2l/3yr y r[1/3] « ~ 3 _ 21 • 

Requiring that f n \ zo ^oo coincides with the experimental value, f n sa 131 MeV, one should take 
a 1 / 3 420 MeV. For such a, the value of 1/a 1 / 3 is close to Zq, i.e., we are in the region oiZq ~ 1 
and we may expect that, even if we use exact formula (|46p with z$ = Zq, the value of f n would 
not change much. Indeed, to get f n sa 131 MeV from Eq. (jloT) for 1/zq = 323 MeV, we should 
take a 1/3 sa 424 MeV = a l J 3 . Thus, in this range of parameters, the value of f n is practically 
in one-to-one correspondence with the value of a. It is convenient to introduce a dimensionless 
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variable 



a = azl = - g^azQ 



(48) 



1 /*3 

Then the values a ' = 424 MeV and = 323 MeV correspond to a = 2.26 = clq. As one can 
see from Fig.([T]), the dependence of f n is practically flat for a > 2. 





FIG. 1: Left: Pion decay constant f n as a function of a for fixed a 1 / 3 = 424 MeV. Right: Function n(a) 



The confinement radius zo presents a natural scale to measure length, so it makes sense to 
rewrite the form factor formula (|38p as an integral over the dimensionless variable £ = z/zq: 



FAQ 2 ) =3 f d(( J(Q,(,z ) 
Jo 



n(a)^ 2 (C,a) + ^^(C,«) 



d(( J(Q,(,zo)p((,a) , 



(49) 



where the mass scale a is reflected by the dimensionless parameter a. The factor n{a) takes care 
of the correct normalization of the form factor. It is given by 

r[2/3] h/M 



ii 



(a) = 2 x /3 a 2/3 



r[l/3] /_ 2/3 (a) 



(50) 



For small a, it may be approximated by |a . For large a, using the fact that l2/z{ a ) / I-2/?>{ a ) 1S 
very close to 1 for a ^ 2, we may approximate n(a) ~ 0.637a 2 / 3 in this region. In terms of n(a), 
the pion decay constant can be written as 



For large a, this gives 



For small a, we have 



\a>2 r 



0.311a 



1/3 



3a 2/3 

— — a 1/3 + . . . ss 0.338 azt, = 0.338 — . 

2y/2lT z 



(51) 



(52) 



(53) 
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The functions (/?(£, a), VKC a) are just the 3> and ^ wave functions written in £ and a variables. 
For a = 0, the limiting forms are (/?(£, 0) = 1 — £ 4 and V'(C)O) = 1- As a increases, both functions 
become more and more narrow (see Fig 12]). For density, we have p{Q,a = 0) = 4£ 2 in the a — ► 




0.2 0.4 0.6 0. 




0.2 0.4 0.6 



FIG. 2: Functions a) (left) and ip(£,a) (right) for several values of a: a = (uppermost lines), a = 1, 
a = 2.26, a = 5, a = 10 (lowermost lines). 

limit, a function that vanishes at the origin (see Fig.([3|)). For nonzero a, the value of p(£ = 0, a) 
monotonically increases with a, and the function itself narrows. The increase of p(( = 0, a) with a 





0.2 0.4 0.6 0.8 



0.2 0.4 0.6 0.S 




0.2 0.4 0.6 



FIG. 3: Le/t: Function p((,a) for a = 0, a = 1, a = 2.26, a = 5, a = 10. Middle: Densities p(£2.26) for 
pion and p p (C) for p-meson in the hard-wall model. Right: Same for densities multiplied by £. 



is generated by the monotonically increasing function n{a). It is interesting to compare the pion 
density p(C, 2.26) (taken at the "experimental" value a = 2.26) with the p-meson density p p (C) °f 



Ref. [23j. These densities are rather close for £ > 0.5, but strongly differ for small £. In particular, 
the p-meson density is more than two times larger for £ = 0, which corresponds to the hard-wall 
model result that g p is essentially larger than f n . 



B. Pion Charge Radius 



It is interesting to investigate how well these values zq = 1/323 MeV and a = (424 MeV) 3 
describe another important low-energy characteristics of the pion - its charge radius. Using the 
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Q 2 -expansion of the vector source 23f] 



j(Q,c,z ) = i- ^zie 



1 - 2 In C 



+ 



and explicit form of the density 



p(C,«) = -r(l/3)r(2/3)a 2 C 4 



Ka)/„ 2/3 « 3 ) 




+ 



/-l/3« 3 

u(a) 



(54) 



z,(a)/ 1/3 (aC 3 ) 



(55) 



where v(a) = \ l2/z{ a )/^-2/ , i{ a )-, we obtain for the pion charge radius: 



{ri) = \4 [dee 



1 - 2 In C 



/3 (C,a) = lz 2 |l-^ + O(a 4 )} 



(56) 




FIG. 4: (r 2 .) in fm 2 for zo = as a function of a. 



Hence, for fixed zq and small a, when a <C the pion radius is basically determined 



by the confinement scale Zq. In particular, (r 2 ) 



| Zq for a 



0. Numerically, taking 



zq = Zq ~ 1/323 MeV = 0.619 fm, we obtain (r 2 ) = 0.51 fm 2 . This result is very close to the value 
(rpk.« 0.53 fm 2 that we obtained in the hard-wall model for the p-meson electric radius determined 



m 



23J from the slope of the Gc{Q 2 ) form factor. However, since Gc(Q 2 ) involves kinematic- type 
terms Q 2 /m 2 , it seems more appropriate to compare F- K (Q 2 ) with the J-\\{Q 2 ) form factor (|42p 



given directly by a wave function overlap integral. The slope of J-u{Q 2 ) is smaller than that of 
Gc(Q 2 ), and the corresponding radius is also smaller: (r 2 )fr = 0.27 fm 2 . Thus, for a = 0, the pion 
r.m.s. radius is about 1.4 times larger than the p-meson size determined by (r 2 )^ 2 . 

With the increase of a, the pion becomes smaller (see FigJH). The experimental value of 0.45 fm 2 
25J is reached for a ~ 0.9. However, the corresponding value f n ~ 80 MeV is too small. If we take 
i = clq = 2.26, then (r 2 ) = 0.34 fm 2 . Thus, if we insist on using zq = Zq dictated by the hard-wall 
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model calculation of the p-meson mass, and the value of a producing the experimental f n (note that 
then a" 2 / 3 « 0.222 fm 2 ), the pion radius is smaller than the experimental value. In linear units, 
the difference, in fact, does not look very drastic: just 0.58 fm instead of 0.66 fm. Given that the 
hard-wall model for confinement is rather crude, the agreement may be considered as encouraging. 
Furthermore, one may expect that, in a more realistic softer model of confinement, the size of the 
pion will be larger. Such an expectation is supported by our soft-wall model calculation of the 
p-meson electric radius, for which we obtained (r 2 p )c =0.66 fm 2 (0.40 fm 2 for (r 2 )jr), i.e., the result 
by 0.13 fm 2 larger than in the hard- wall model. If (r 2 ) would increase by a similar amount, the 
result will be very close to the quoted experimental value. 

To find (r 2 ) for large a (i.e., when a •> z^~ 3 for fixed zq, or when zq > a" 1 / 3 for fixed a), we use 
first the observation that, in the region a > 2, we may approximate v(a) ~ 1. Then the factor in 
square brackets in Eq. (f55l) becomes a function of the combination aC 3 = /i (call it and we 

can write 

2/3 



(4) 



^ 3 

o>2~ 4 



r(l/3)r(2/3) (-) J d/V /3 i2(A0 



1 



3 a 



(57) 



For a ^ 2, the upper limit of integration in this expression may be safely substituted by infinity 
producing 

2 2/3 



3T 2 (2/3) 



= G 



d/i n 5/3 R(ji) In n w G In 0.566 



(58) 



which gives 



(rl) 



r(i/3) 



^2/3 



1 + - In , 

3 V0.566 



a>2 2 4 /3r(2/3) \aj 

Using Eq. (|47p . we can express the coefficient in front of the square bracket in terms of /„■: 

3 1 



K) 



+ 



In 



az 



o 



(59) 



(60) 



a>2 4vr 2 / 2 2vr 2 / 2 V 0.566, 
Thus, (r 2 ) in the a ^ 2 region consists of two componens: a fixed term 3/47T 2 / 2 and a term 
logarithmically increasing with zq. As zq — > oo, the pion charge radius becomes infinite, reflecting 
the fact that the pion in this model is massless. A similar structure in the expression for the pion 
charge radius was obtained 



281 ] in the Nambu-Jona-Lasinio (NJL) model 



+ 



1 



In 



mt 



(61) 



It also has the logarithmic term 
the infrared-finite piece 3/27T 2 / 2 



27T 2 / 2 8^ 2 / 2 

nm 2 resulting in the infinite radius for massless pion and 
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32l | . The latter, however, is twice larger than that in our result 
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(|60p and contributes 0.34 fm 2 to (r 2 ), with the chiral logarithm term producing the extra 0.11 fm 2 
required for agreement with experiment. In our case, the logarithmic term taken for a = ag is 
approximately equal to 3/4.7T 2 / 2 , thus almost doubling the outcome value for (r 2 ). More precisely, 
we can write 

3 



1 , ( a 

1 + - In 

3 V2.54 



(62) 



a>2 27T 2 / 2 

For a = 2.26, the modified logarithmic term gives a very small contribution, and our net result is 
very close to the value given by the NJL fixed term. Numerically, though, this prediction of the 
hard- wall AdS/QCD model, as we have seen, is essentially smaller than the experimental value. 



C. Form Factor at Large Q 2 



In the large-Q 2 limit, the source J7"(Q, z) is given by its free-field version zQK\(Qz) that behaves 
asymptotically like e~® z . As a result, only small values z ~ 1/Q are important in the form factor 
integral, and the large-Q^ asymptotic behavior of the form factor is determined by the value of 



p(z) at the origin 



24 ]. namely, 



2p(0) 2^(0) 4tt 2 / 2 



F n (Q 2 



Q 2 



(63) 



Note that the combination 47r 2 / 2 



so 



Q 2 Q 2 Q 2 

0.68 GeV 2 frequently appears in the pion studies. In 



particular, it is the basic scale of the pion wave function in the local quark-hadron duality model 
33, 



341 ] , where it corresponds to the "pion duality interval" . 




Q 2 F n (Q 2 ) 




10 



15 



20 



Q 



FIG. 5: Left: Contributions to pion form factor F 7r (Q 2 ) from \l/ 2 -term (lower curve), from <I> 2 -term (middle 
curve) and total contribution (upper curve). Right: Same for Q 2 F 7r {Q 2 ) . 



The leading contribution comes entirely from the $ 2 term of the form factor integral (|41|) while 
the ^ 2 term contribution behaves asymptotically like 1/Q 4 since it is accompanied by extra z 2 
factor. Note, however, that it is quite visible in the experimentally interesting region Q 2 < 10 GeV 2 : 
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it is responsible for more than 20% of the form factor value in this region (moreover, at Q 2 = 0, 
the ^ 2 term contributes about 40% into the normalization of the form factor). 

From a phenomenological point of view, different AdS/QCD-like models for the pion form factor 
differ in the shape of the density p(() that they produce. If we require that the density p(z) equals 
2ir 2 f 2 at the origin, the asymptotic behavior is F n (Q 2 ) — > sq/Q 2 in any such model. For Q 2 = 0, 
the form factor is normalized to one, so basically the models would differ in how they interpolate 
between these two limits. In particular, the simplest interpolation is provided by the monopole 
formula 



\Q 2 ) 



i 



(64) 



i + Q 2 /s ' 

while our hard- wall calculation gives a curve that goes above F™ ono (Q 2 ): the ratio 
F 7r (Q 2 )/F™ no (Q 2 ) is larger than 1 for all Q 2 > 0, slowly approaching unity as Q 2 — ► oo (see 

FigED. 



FAQ 2 ) 




FAQ)IF„ mon °(Q) 
1.2 




20 40 60 80 100* 



FIG. 6: Left: Pion form factor F V (Q 2 ) from the holographic model (upper curve) in comparison with the 
monopole interpolation F™ no (Q 2 ) (lower curve). Right: Ratio F^{Q 2 )/F^ ono (Q 2 ). 



In fact, a purely monopole form factor was obtained in our paper [24j, where we studied the 
/9-meson form factors in the soft-wall holog rap hic model, in which confinement is generated by ~z 2 



(65) 



oscillator-type potential. It was shown in [24 ] that the form factor integral 



HQ 2 ,*) 



dzzJ°(Q,z)\<5>(z,K) 



in which <3?(z, k) = ^/2ne z2k2 / 2 is the lowest bound state wave function, and 



J°(Q,z)=z 2 . 2 f- 
Jo U 



dx 



1 — X 



2 2 

Z K 



(66) 



is the bulk-to-boundary propagator of this oscillator-type model, is exactly equal to 
1/(1 + Q 2 /(4k 2 )). The magnitude of the oscillator scale k was fixed in our paper [24] by the 
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value of the p-meson mass: k = k p = m p /2. As a result, the form factor ^{Q 2 , k = m p /2) had the 
p-dominance behavior 1/(1 + Q 2 /m 2 ). 

If we take k = K n = irf n « 410 MeV both for &(z,k) and J°(Q,z), the integral (f65|) 
gives 1/(1 + Q 2 /s ). The relevant wave function &(z,K n ) has the expected correct normaliza- 
tion <E>(0, Kjt) = V2-irf w , however, the slope 1/sq of 1/(1 + Q 2 /sq) at Q 2 = (corresponding to 
0.35 fm 2 for the radius squared) is smaller than that of the experimental pion form factor. Fur- 
thermore, Q 2 F™ ono (Q 2 ) tends to s ps 0.68 GeV 2 for large Q 2 , achieving values about 0.5 GeV 2 



26J mea- 



for Q 2 ~ 2 GeV 2 , and thus exceeding by more than 25% the experimental JLab values 
sured for Q 2 = 1.6 and 2.45 GeV 2 . The authors of Ref. [^tJ proposed to use Eqs. (J65l) . (f66j) as an 
AdS/QCD model for the pion form factor, with k = 375 MeV chosen so as to fit these high-Q 2 
data. However, such a choice underestimates the value of f 2 by almost 30%. Our opinion is that 
the AdS/QCD models should describe first the low-energy properties of hadrons, and the basic 
low-energy characteristics, such as m p and /„-, should be used to fix the model parameters. On the 
other hand, if the form factor calculations based on these parameters disagree with the large-Q 2 
data, it is quite possible that this is just an indication that one is using the model beyond its 
applicability limits. Furthermore, as we have seen in the hard-wall model, to correctly describe 
the pion one needs to include the chiral symmetry breaking effects absent in the vector channel. 
As a result, equations for pion wave functions are rather different from those in the p-meson case. 
Similarly, there are no reasons to expect that, in a soft- wall model, the pion density should have the 
same shape as the p-meson one. Unfortunately, the procedure of bringing in the chiral symmetry 



breaking effects that was used in the hard-wall model of Ref. 



(J faces serious difficulties when 



) model with the z 2 soft wall. As discussed in Ref. ji]], the solution of 
the equation for the X field in this model requires that chiral condensate a and the mass parameter 
m q are proportional to each other, so that a cannot be varied independently of m q . Moreover, if 
one takes the chiral limit m q = 0, the chiral condensate should also vanish. This difficulty may 
be avoided by switching to more sophisticated recent models (cf. 16l. 1 171. Ilq) in which the chiral 
condensate is generated dynamically. However, such a consideration goes well beyond the scope of 
the present paper. Thus, we just resort to an idea that whatever the mechanism is involved, the 
net practical outcome is a particular shape of the density p(z) that eventually determines the pion 
form factor and other pion characteristics. Below, we give an example of a density p mod (z) that 
is normalized at the origin by the experimental value of /„-, i.e., p mod (0) = 2k 2 , but which is also 
capable to reproduce the experimental value of the pion charge radius. 

Evidently, to increase the radius, we should take a density which is larger for large z than 
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7T /„■). Since the overall integral normalization of the density is kept fixed, this can be 



achieved only by decreasing the density for small z values. 



z p(z, a) 
1.75 




mod (e 2 )/Fj mo n o (e 2 } 



0.98 

0.96 
0.94 
0.92 
0.9 
0.88 



10 



15 



20 



Q 1 



FIG. 7: Left: Model density zp mod (z) (measured in fm 1 ) is larger than the density z|$(z, n-n)\ 2 for large z 
(displayed in fm). Right: Ratio F™ od (Q 2 )/F™ ono (Q 2 ) for B = 1/4. 



Consider a simple ansatz (see FigfTJ) 



„mod 



(z) = 2k\ e"* 2 "? [1-A z 2 K 2 p + B z 4 k 



(67) 



with A = 1 — k 2 /k 2 + 2B. It has both the desired value for z = and satisfies the normahzation 



condition 



dzzp mod (z) = 1 



(68) 



Integrating it with J (Q, z) taken at k = n p produces the model form factor given by the following 
sum of contributions of the three lowest vector states: 



r>mod 



(<5 2 



2 - (1 - 2B)s /m 2 l-(l-4B)s /m 



1 + Q 2 /m? p 



1 + Q 2 /2mj 



+ 



2BsQ/m 



1 + Q 2 /3m 2 



The slope of F™ od (Q 2 ) at Q 2 = is given by 



dF™ d {Q 
dQ 2 



1 



(69) 



(70) 



Taking B = 1/4, one obtains the experimental value 0.45 fm 2 for (r 2 ). It is interesting to note 
that the model density providing this value, has an enhancement for larger values of z (see FigfT]), 
just like the pion densities in the hard- wall model (see Fig.Q). Due to a larger slope, F™ d (Q 2 ) 
decreases faster than the simple monopole interpolation F™ ono (Q 2 ) and, as a result, is in better 
agreement with the data. In fact, it goes very close to Q 2 < 1 GeV 2 data, but exceeds the values 
of the JLab Q 2 =1.6 and 2.45 GeV 2 points by roughly 10% and 20%, respectively. 

This discrepancy has a general reason. The asymptotic AdS/QCD prediction is 
Q 2 F n (Q 2 )\Q2^ oa — > 47T 2 / 2 which is « 0.68 GeV 2 for experimental value of f n . On the other hand, 
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JLab experimental points correspond to Q 2 F^ xp (Q 2 ) ~ 0.4 GeV 2 , which is much smaller than the 
theoretical value quoted above. The pre-asymptotic effects, as we have seen, reduce the discrep- 
ancy, but there still remains a sizable gap. As we already stated, such a disagreement may be just 
a signal that we are reaching a region where AdS/QCD models should not be expected to work. In 
particular, AdS/QCD models of Refs. P, U, If], s| describe the pion in terms of an effective field or 
current, without specifying whether the current is built from spin- 1/2 fields, or from scalar fields, 
etc. For Q 2 above 1 GeV 2 , the quark substructure of the pion may be resolved by the electromag- 
netic probe (which is a wide-spread belief), and the description of the pion "as a whole" may be 
insufficient. 

V. SUMMARY 

In this paper, we studied the pion in the chiral limit of two flavor QCD. To this end, we described 
a formalism that allows to extract pion form factor within the framework of the holographic dual 
model of QCD with hard- wall cutoff. Following Ref. Q|, we identified the pion with the longitudi- 
nal component of the axial-vector gauge field. We defined two (Sturm-Liouville) conjugate wave 
functions &(z) and ^f(z) that describe the structure of the pion along the 5 th dimension coordinate 
z. These wave functions provide a very convenient framework to study the holographic physics of 
the pion. We demonstrated that, just like in the p-meson case 23[], the pion form factor is given 



by an integral involving the function p(z) that has the meaning of the charge density inside the 
pion. However, in distinction to the p-case, when the density was simply given by \$(z)\ 2 , the pion 
density has an additional term proportional to ^(z)) 2 and entering with the z-dependent coeffi- 
cient reflecting the mechanism of the spontaneous symmetry breaking. Both terms are required 
for normalization of the form factor at Q 2 = 0. 

We found an analytic expression for the pion decay constant in terms of two parameters of the 



model: a and zq, similar to those used in Ref. 



. Analyzing the results, we found it convenient 



to work with two combinations a = g$o /3 and a = <xzq of the basic parameters. In particular, we 
found a = ao = 2.26 for the value of a corresponding to the experimental p-meson mass m p and 
pion decay constant /„-. The importance of the parameter a is that its magnitude determines the 
regions, where the pion properties are either governed by the confinement effects or by the effects 
from the spontaneous chiral symmetry breaking. For example, in the practically important domain 
a > 2, the pion decay constant is determined primarily by a, with negligibly tiny corrections due to 
zq value. However, when a < 1 the pion decay constant is proportional to the ratio a/ zq. Besides, 
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for small a <C 1, the radius of the pion is given by (r 2 ) = |zq, i.e., as one may expect, the pion 
size is completely determined by the confinement radius. On the other hand, for a > 2 the radius 
is basically determined by l/c 1 / 3 , slowly increasing with zq due to the lna/ao correction. 

We also found that the pion rms charge radius (r 2 ) 1 / 2 ~ 0.58 fm in the hard- wall model is smaller 
than that measured experimentally. In a sense, the hard wall at the distance zq ~ 0.62 fm (fixed 
from the p- meson mass), "does not allow" the pion to get larger. So, we argued that if the IR wall 
is "softened" , the size of the pion may be increased by an amount sufficient to accomodate the data. 
A straightforward idea is to use the soft-wall model of Ref. and treat the pion in a way similar 
to what was done in [24] for the p- meson case. Unfortunately, there are prohibiting complications 
with directly introducing the chiral symmetry effects within the AdS/QCD model with the z 2 soft 
wall. As explained in Ref. {|, the chiral condensate a in such a model is proportional to the mass 
parameter m q , so that in the chiral limit the condensate vanishes together with the quark mass. 

To illustrate a possible change in the form factor predictions due to the softening of the IR wall, 
we proposed an ansatz for the pion density function and used the vector current source from the 
soft- wall model considered in Ref. 24|. We demonstrated that this ansatz is capable to fit the 
experimental value of the pion charge radius. It also closely follows the data in the Q 2 < 1 GeV 2 
region, while still overshoots available data in the Q 2 ~ 2 GeV 2 region. The basic source of 
this discrepancy is very general: the asymptotic AdS/QCD prediction for the pion form factor 
is Q 2 F 7V (Q 2 ) — ► 4-7T 2 / 2 , and if one takes the experimental value for /„., one obtains Q 2 F n (Q 2 ) — ► 
0.68 GeV 2 , which is much larger than the 0.4 GeV 2 value given by Q 2 ~ 2 GeV 2 JLab data. For 
this reason, we argued that the disagreement mentioned above may be a signal that the region 
Q 2 > 2 GeV 2 is beyond the applicability region of AdS/QCD models. 



351 ] addressing the same 



Finishing the write-up of this paper, we have learned that the paper 
problem was posted into the arxive. We did not observe, however, essential overlaps with our ideas 
and results. 
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